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^ Abstract 

We consider the problem of computing a response curve for binary cellular automata - that 
is, the curve describing the dependence of the density of ones after many iterations of the rule 
on the initial density of ones. We demonstrate how this problem could be approached using rule 
130 as an example. For this rule, preimage sets of finite strings exhibit recognizable patterns, 
and it is therefore possible to compute both cardinalities of preimages of certain finite strings 
and probabilities of occurrence of these strings in a configuration obtained by iterating a random 
J> initial configuration n times. Response curves can be rigorously calculated in both one- and 

C — two-dimensional versions of CA rule 130. We also discuss a special case of totally disordered 

initial configurations, that is, random configurations where the density of ones and zeros are 

^ equal to 1/2. 

00 

O Introduction 

Cellular automata (CA) can be viewed as computing devices, which take as an input some initial 
configuration. The CA rule is iterated a number of times starting from this configuration, resulting 
in a final configuration, which constitutes the output of the computation. If the CA rule is complex, 
the above computation may be very difficult to characterize, let alone understand in detail. If one 
considers "simple" CA rules, however, one can say quite a lot about the process, as we shall see in 
what follows. 

In many practical problems, e.g., in mathematical modelling, one wants to know how a CA rule 
iterated over an initial configuration affects certain aggregate properties of the configuration, such 
as, for example, the density of ones. If we take a randomly generated initial configuration with a 
given density of ones, and iterate a given rule n times over this configuration, what is the density 
of ones in the resulting configuration? Using signal processing terminology, we want to know the 
"response curve", density of the output as a function of the density of the input. Response curves 
appear in computational problems, and a classical example of such a problem in CA theory is the 
so-called density classification problem (DCP). If we denote the density of ones in the configuration 
at time n by c n , the DCP asks us to find a rule for which Coo = 1 if Co > 1/2 and Coo = if cq < 1/2 
- that is, CA rule with density response curve in a form of a step function. Since it is known that 
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such a rule does not exist (Land & Belew, 1995), once could ask a related general question: which 
response curves are possible in CA rules? Obviously, this problem is much more difficult than DCP, 
and not much is known about it. We propose to approach this problem from an opposite direction: 
given the CA rule, what can we say about its response curve? It turns out that in surprisingly 
many cases, the response curve can be calculated exactly, providing that preimage sets of finite 
strings under the CA rule exhibit recognizable patters. We will demonstrate this technique using 
as an example elementary CA rule 130, in both one- and two-dimensional spaces. 



Basic definitions 

Let M = {0, 1} be called a symbol set, and let S(Af) be the set of all bisequences over M , where by a 
bisequence we mean a function on % to N. Throughout the remainder of this text the configuration 
space S(N) = {0, 1} Z will be simply denoted by S. 

A block of length n is an ordered set b$b\ . . . 6„_i, where n £ N, b\ 6 M . Let n 6 N and let B n 
denote the set of all blocks of length n over Af and B be the set of all finite blocks over J\f. 

For r € N, a mapping / : {0, l} 2r+1 i->- {0, 1} will be called a cellular automaton rule of radius r. 
Alternatively, the function / can be considered as a mapping of ^2r+i into Bq = N = {0, 1}. 

Corresponding to / (also called a local mapping) we define a global mapping F : S — > S such 
that (F(s))i = f(si- r , . . . , Si, . . . , Si +r ) for any s £ S. The composition of two rules f, g can be now 
defined in terms of their corresponding global mappings F and G as (F o G)(s) = F(G(s)), where 
s e S. 

A block evolution operator corresponding to / is a mapping f : B i-} B defined as follows. Let 
r 6 N be the radius of /, and let a = a^ai . . . a n -i £ B n where n > 2r + 1 > 0. Then 

f (a) = {f(a t , a l+1 , a i+2r )}^ r -\ (1) 

Note that if b € B 2r +i then f(b) = f (b). 

We will consider the case of N = {0, 1} and r = 1 rules, i.e., elementary cellular automata. In 
this case, when b G B 3 , then f(b) = f(6). The set B 3 = {000,001,010,011, 100, 101, 110, 111} will 
be called the set of basic blocks. 

The number of n-step preimages of the block b under the rule / is defined as the number of 
elements of the set f _n (6). Given an elementary rule /, we will be especially interested in the 
number of n-step preimages of basic blocks under the rule /. 

As mentioned in the introduction, we will use as an example rule 130 (using Wolfram's num- 
bering scheme (Wolfram, 1994)) with local function defined as 

f( rvz \-l 1 if (zj/z) = (001)or(lll), (2) 
n x 'y> z )-\ otherwise. [Z) 

From now on, both the local function / and the corresponding block evolution operator f will refer 
to rule 130, unless otherwise noted. We will calculate the response curve for this rule by considering 
the structure of preimages of finite strings, using the approach described in earlier papers (Fuks, 
2006; Fuks & Haroutunian, 2009; Fuks, 2010). 



Structure of preimage sets in one dimension 

Since the only preimages of 1 under rule 130 are blocks 111 and 001, let us consider first their 
preimages, that is, the sets f — "(111) and f _ri (001). The following two propositions describe these 
sets. 
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Proposition 1. The setf "(111) has only one element, namely the block 11 . . . 1 , hence cardf "(111) = 

2n+3 

1. 

Proof. Exhaustively checking all 32 potential preimages of 111 one can show that the only length-5 
string b such that {(b) = 111 is 11111. Since any block of l's consists of overlapping blocks 111, its 
preimage also must consist entirely of ones, and the proof by induction follows. □ 

Proposition 2. The set f~ n (001) consists of all blocks of the form 

10 1 (if i is odd) or <L^l£ 1 (if i is even), 

2n-2i 2i 2n-2i 2i 

where i € {0 . . . n} and * denotes an arbitrary value in S. 

Proof. By induction, if n = 1, then i £ {0, 1}, and our formula provides us with the following two 
types of preimages, **001 or 10111. One can exhaustively check all 32 blocks of length five to 
verify that these are the only five blocks b such that {(b) = 001. 
Now, let us assume that we have the following n-step preimage: 

1 where % even and i G {2 . . . n}. 

2n-2i 2i 

We find the preimages of this string, starting from the right and working toward the left. By 
considering all 16 blocks of length four, we can see that the only string that has image 11 under f 
is 1111, so we start with the following preimage (written above the string we are considering, with 
arrow indicating direction of proceeding) : 

<- 1 1 1 1 
**...** 1 1 ... 1 1 1 

Continuing to the left, using the rule table of rule 130 we can construct the preimage up to the 
following point: 

<- 1 1 1 1 ... 1 1 1 1 
*★...** 1 1 ... 1 1 1 

From here to the left, all the remaining entries in the preimage are arbitrary: 

★ ★ ... * * 1 1 1 1 ... 1 1 1 1 
★ * ... * * 1 1 ... 1 1 1 

We now have the following (n + l)-step preimage: 

★ 101 where (i + 1) is odd and (i + 1) G {3 . . . (n + 1)}. 

2(n+l)-2(i+l) 2(i+l) 

To finish the proof, we need to perform similar analysis for two other cases as follows. 

• When i = 0, a similar argument is used to find all (n + l)-step preimages in which i = 0, 1. 

• When i € {1, . . . n} and odd, another similar argument can be used to find all (n + l)-step 
preimages in which (i + 1) is even. 
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We omit the details, but note that all possible (n + l)-step preimages are accounted for in the 
analysis, thus completing the induction step. □ 

We may now proceed with enumerating the elements of /~ n (001). 
Proposition 3. For rule 130 in one dimension, we have cardf -n (001) = . 



Proof. According to Proposition (2), for each value of i £ {0. . . n}, there are 2n — 2i arbitrary 
values in the preimage, so that there are 2 2n ~ 2t of such preimages. Summing over i we obtain 

n n 4 n+l _ 1 

card[/-"(001)] = 2 2 "~ 2i = 4 ™ Yl 4 ~* = § ■ ( 3 ) 

i=0 i=0 

□ 

Similar reasoning as in Proposition (1) leads to 

cardr n (110) = cardf- n (101) = 1. (4) 

It turns out that we now have enough information to find the cardinality of the preimage sets of 
each of the four remaining basic blocks, since not all of them are independent. The interdependence 
between preimage sets can be easily understood if one considers this problem using the language of 
probability theory. Suppose that we start with an infinite configuration in which all sites indepen- 
dently assume the value or 1 with equal probability 1/2. If we iterate our rule n times starting 
with this initial condition, one can show (Fuks, 2010) that the probability P n (b) of occurrence of 
block b in the final configuration is given by 

P n {b) =2-l 6 l- 2n cardf- n (6), (5) 

where \b\ is the length of the block b. Probabilities of different blocks are not independent, since 
the following consistency conditions (Dynkin, 1969) hold: 

P„(000) + P n (001) = P„(100) + P n (000) = P n (00), (6) 

Pn(H0) + P n (lll) = Pn(OH) + P„(1H) = P„(ll), (7) 

P n (010) + P n (110) = P n (100) + P n (101) = P n (10). (8) 
Equation (6) implies that P„(100) = P„(001), hence 

cardf- n (100) = cardf- n (001). (9) 
Likewise, equations (7) and (8) imply that 

cardf- n (011) = cardf- n (110), (10) 

cardf- n (010) = cardf -n (101) + cardf" n (100) - cardf _n (110). (11) 

Finally, since we know the total number of all preimages of all basic blocks, we may subtract the 
seven known formulae from the total to find card f~ n (000). The results are summarized in Table 1. 
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b 


cardf _A/ (6) 


b 


card[/ n (b)} 





4 n+l _ 3 


1 


3 


1 


3 


1 1 


1 


1 


3 


1 1 


1 


1 1 


1 


1 1 1 


1 



Table 1: Number of Preimages of Basic Blocks of One-Dimensional Rule 130 



Dependence on the initial density in one dimension 

Let us now define p to be the probability that a given cell is in state 1 in the randomly generated 
initial configuration. Obviously, 1 — p is then the probability that a cell is in state 0. Since in the 
initial configuration all cells are independent, the probability that a given block b will occur in the 
initial configuration is 

P„(6) = p# of rs in \l - P )# of °' s in b . (12) 

Thus, if we wish to find the density of ones after n iterations of rule /, we can use preimages of 1 
and write 

P n (l)= P o(6). (13) 

6e/-»(i) 

Detailed derivation and discussion of the above equation can be found in (Fuks, 2010), and will not 
be repeated here. Now, since the only preimages of 1 are 001 and 111, one obtains 

p n (i) = p n _i(ni)+p„_i(ooi)= p °( fe )+ E p o^- ( 14 ) 

6gf-»+l(ni) 6Gf- n+1 (001) 

We know the structure of the n-step preimage sets for Rule 130 in one dimension, thus we may 
compute desired probabilities in the above formula. We know from Proposition 1 that the preimage 
of 111 is comprised entirely of ones, hence 

p„(iii) = P (n^_i) = P 2n+3 . (15) 

2n+3 

From Proposition 2 we find, 

p n (ooi)= (i-p)V m + E ( i -p)p 2i+2 (is) 

ie{0...n},i is even ie{0...n},i is odd 

= E a - P ) 2 P 4k+1 + E a - p)p Ak+A (17) 

k=0 k=0 

p (_p4f(n-l)/2l+4 + p 4r(n-l)/2l+5 _ p 4L(n+l)/2j+3 + p 3 _ p + ^ 



p 3 + p 2 + p+l 



(18) 



Using the fact that P n (l) = Pi-i(Hl) + Pt-i(001), we obtain the equation of the response curve, 
that is, dependence of the density of ones after n steps, denoted by c n = P n (l), on the initial 
density, denoted by p: 

p (_ p 4r(n-2)/2l+4 + p 4r(n-2)/2l+5 _ p 4Ln/2j+3 + p 3 _ p + ft 

C n = p n+ + 3|2l Tl • ( 19 ) 

P 6 + p Z + p + 1 
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Figure 1: Plots of the response curve, i.e., dependence of Cqo on p for one-dimensional rule 130. 
Circles correspond to computer simulations using a lattice of 5000 sites, as described in the text, 
while the continuous line represents theoretical curve. 



By taking the limit of n 



oo, we obtain the asymptotic response curve, 
p{p 3 -p + l) 



p 3 + p 2 + p + 1 
1 if p 



if p< 1, 
1. 



(20) 



We performed computer simulations to illustrate the formula (20) for the the dependence of Cqq 
on the initial density p. We considered an initial configuration of 5000 cells, and varied the initial 
density from to 100%, increasing it by a step size of 1%, iterating rule 130 until we reached a fixed 
density. Results were then averaged over 20 runs for each initial density. The results are presented 
in Figure 1. As we can see, the response curve calculated above for infinite configurations agrees 
very well with simulations performed on a finite lattice. As a final remark in this section, let us 
note that for the special case of p = 1/2 we obtain 



1 + 14-. 
6 3 



(21) 



We can see that the convergence toward Cqo = 1/6 is exponential. 

We will now proceed to consider two-dimensional version of rule 130. As it turns out, the basic 
ideas used in the preceding section for calculating response curve in one dimension can be carried 
over to two dimensions. 



Basic definitions in two dimensions 

All basic definitions introduced at the beginning of this paper can be easily generalized to two 
dimensions. We will consider a variant of rule 130 in which the state of the cell at the next 
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iteration depends on the cell itself, its right neighbour and the top neighbour. We will call this 
type of neighbourhood an L-shaped neighbourhood. To be more precise, let us define the local 
mapping of a two-dimensional cellular automaton with L-shaped neighbourhood as 

g : (* z ) -+ {0,1}, (22) 

where x,y, z G G. Similarly as in one dimension, g has a corresponding global mapping, G : 
{0, l}^ 2 -»• {0, l}^ 2 such that 

(G(s)). j = 9 ( s ^ Si+h3 ) (23) 

for any s G {0, l}^ 2 . Blocks in two dimensions will be defined as regions of 2D lattice in the shape 
of isosceles right triangles. We define the set of triangular blocks of size r as the set consisting of 
elements 



b 



(24) 



bi,i . . . b r ,i, 

where each bij G {0, 1}. This set will be denoted by %■ The block evolution operator g : % — > T r -i 
will be defined as a function which transforms triangular block (24) into another block 



(25) 



Cl,r- 1 
Cl,l . . . Cr-1,1, 

where c i;j = <l{'' h \' b . +lj ) for i G {1, . . . , r - 1}, j G {1, . . . , r - i}. 
Finally, rule 130 is defined in two dimensions as 

Kv)-\ otherwise. 1 bj 

Since the only preimages of 1 are triangular blocks (q i) and (} 1 ), we need, similarly as in one 
dimension, to analyze structure of preimage sets g~™( [j i ) and g~ n ( \ 1 ). 

Structure of preimage sets in two dimensions 

Proposition 4. The set g _n ([j i ) consists of all blocks of the following form: 



ai 

a 2 1 

a 3 : ■ ■ 

-k ... ★ aj+2 1 1 



(27) 



i+l 
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where i 6 {0 . . . n} and for 03 . . . 0^+2 we take arbitrary values. Ifi = 0, then (a±, 02) = (0, 0), w/ii/e 
ifi>0, 01,02 are determined by 



i+l 

I 



1 + ( ■ ! K' mod2 > ( 28 ) 



ai 

i+l ✓ 

a 2 = 1 + 5^ I * )aj + i mod2, (29) 

Before we attempt the proof, we need to observe two facts. First of all, note that 

g(yi)=x + y+l mod 2, (30) 

and let us define h(x, y) = x + y + 1 mod 2. Secondly, consider the following procedure. We start 
with a binary sequence a±, 02, ...Ofc and replace it by a sequence of pairs h(ai, 02), ^(02, 03), . . . , h(a k -i, a k ). 
This new sequence is obviously of length A; — 1. If we repeat this process k — 1 times, we will end 
up with just one number, to be denoted by q(ai, 02, ■ ■ ■ , a k ). The following lemma can be easily 
proved by induction. 

^ ( k l\ 

Lemma 1. q(ai, 02, . .., a k ) = 2 k ~ 1 — 1 + ^ ( . -|) a * m °d 2. 

Proof. When k = 2, one obviously has g(ai, 02) = h(a\, 02) = «i + 0,2 + 1 mod 2, and 

- 1 + Yl (■ 1 i) a * mod2 = l + ai + a 2 mod 2, (31) 

thus the lemma is indeed true for k = 2. 

Suppose now that it holds for a given k, and let us consider a binary sequence a±, 02, . . . , o& + i. 
We clearly have 

g(ai, a 2 , . • • , a fc+ i) = g(ai, a 2 , . . . , a k ) + q(a 2 , a 2 , . . . , a fc+ i) + 1 mod 2, (32) 

hence 

9(01,03, . . . , a k+1 ) = 2 k ~ l - 1 + ( k Z I) a * + " 1 + E 1 2) °* + 1 m ° d 2 (33) 

= 2^ - 1 + ai + ( • r 1 ) ai + afe +! + e ( ■ 1 2 ) ai mod 2 ' 



(34) 



i=2 v ' i=2 

Using Pascal's identity, this becomes 

q(ai,a 2 , ... , a k +i) = 2 k - 1 + ai + a k + ^ f _ j mod 2. 
Using (q) = (fc) = 1 we can incorporate ai and afc + i into the sum, obtaining 



(35) 



k+l ( k \ 

q(ai,a 2 , a k+1 ) = 2 fc -l + ^(._ 1 j mod 2 ' 

i=i ^ ' 



which is the desired formula for At + 1. Proof by induction is therefore complete. 



(36) 



□ 
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Having this lemma, we can now proceed with the sketch of the proof of the Proposition 4. 
Proof. Consider first the case of n = 1, when the only preimages we obtain are of the form 

★ cti 

★ or a 2 1 (37) 

★ 1 3 1 1 

where a\ = 03 and a 2 = 1 + 03 mod 2 are given by eq. (28) and (29) for any 03 € {0, 1}. One can 
check that these are indeed the only desired preimages by applying g to all blocks T3 and verifying 
that only the blocks of the above form produce [j 1 . This means that the block [j 1 can appear in 
a configuration by two ways, by moving from the left one unit at a time, or by being created from 
a configuration shown on the right hand side of (37). 

Consider a number of steps larger than one, denoted by n. A given block may appear in a 
certain place because it was created in another location, k units to the right, and then moved to 
the desired place from there in k steps, where k < n, or because it appeared in this place as a result 
of the second configuration of eq. (37). 

Now note that the block shown in eq. (27) in i iterations will produce block 



(38) 



: 

? 1, 



where ? denotes some value resulting from iteration of the rule. This is because the triangular 
block of ones in eq. (27) will shrink until it becomes a single 1, as shown above, and the column 
ai, 02, . . . a m will also shrink until it becomes a pair (u, v), where u = q(ai,a 2 , ■ ■ ■ , Oj+i) and v = 
q(a 2 , as, ... , ai +2 ). The reason for this is the fact that iteration of the rule is equivalent to iterative 
application of h(x,y) to 01,02, . . . Oj+2- Conditions (28) and (29) are equivalent to 



i+l 



1 = E(,-_iH- mod2 ' ( 39 ) 
j= i \ j / 

i+l , . X 

1 = E( ■ * mod2 ' ( 4 °) 



and by Lemma 1 this implies u = 0, v = 0, that is, we obtain column of two zeros, as shown in eq. 
(38). 

Iterating block shown in eq. (38) further, after n — i iterations, § 1 will move to the left one 
step at a time. In the end, as a result of applying operator g n-times to block (27), [j 1 will be 
produced. Since i varies from to n this means that indeed all possibilities of arriving from some 
place on the right (i = 1, . . . n — 1) as well as being "created in place" (i = n) are covered. □ 

Using very similar argument, one can prove analogous proposition for block \ 1 . 
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Proposition 5. The set, g n (}j), consists of all blocks in the set \J™ =0 (Ai \ Bi) U C, where, for 
a fixed value of i, Ai and Bi are, respectively, the sets of all blocks of the form 



bi 

b 2 1 
b 3 '■ 

* b i+ 2 1 ■ 



i+l 



and C is the set whose only element is the block 



: * 

: : 1 

★ . . . * 1 ... 1, 

n-i i+2 



n+2 



where i G {0 . . . n} and for 63 . . . 6^+2 we take arbitrary values. If i = 0, then (b\, 62) = (1, 1) and if 
i > 0, then b\, 62 are determined by 



i+l 



3=2 
i+l 



fo i = E(/ 1 ) h i mod2 > 



b 2 = Yj\j -ip +i mod2 - 



(41) 



(42) 



Proposition 6. The number of ' n- step preimages of blocks [j 1 and \ 1 is given by 

2 n 

card[g-(0 i )]=2^^2-^ i (43) 

i=0 

«Hg-»( i , )] = 2^ (4 - 3 x; 2-^) , (44) 

i=0 

Proo/. For a given value of i, we have (n + 2) + (n + 1) H h(i + 3)+« = w(w+5) 2 " t( ' +5) + i arbitrary 



n(n + 5) — i(i + 5) 



values represented by * in the preimage of § 1 . Therefore, for each i, there are 2 2 

re(rt+5) — r(z+3) 

2" 2 - possible configurations. Summing over all values of i G |0...n} gives the desired 
result. 
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n(n+5)-i(i+3) n(n + 5)-i(i+5) 

In preimages of } 1 , for a given i we have 2 2 ~~ possible blocks in set Ai, and 2 2 
blocks in set B{. Summing over i and adding one block consisting of all ones, we obtain 

n n 

,r -n/1 \1 \ ^ n(n + 5)-a(»+3) x ^ n(n + 5)-»(i + 5) 

card[g "({J] =JJ2 2 _J^2 2 + 1. (45) 

i=0 i=0 



After changing the summation index i in the second sum to j = i + 1 , both sums can be combined 
together, and after simplification one obtains 

2 n 

card[g-«({ 1 )]=2^( 4 -3j:2-^). (46) 

i=0 

□ 

Similar analysis can be performed for block \ , leading to the same number of preimages as in 
the case of \ x , 

2 n 

card[g-«(} )] =2^(4-3^2-^). (47) 

i=0 

There are two other blocks which have easy to describe preimage sets, as demonstrated in the next 
two propositions. 

Proposition 7. The set g~ n ( 5 1 ) consists of all blocks of the following form: 



ai 

a 2 1 

as : '•• ( 48 ) 

a n +2 1 1 



n+l 



where for a% . . . a n+ 2 we take arbitrary values. Ifn = 0, then (ai, 02) = (0, 1), while if n > 0, a±, 02 
are determined by 



n+l 

n 



3=2 

n+l 



ai = i + ^^j iJ a J mod2 ' ( 49 ) 



0,2 = ( . A a j+i mod 2. (50) 

j=2 ~ / 
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Proposition 8. The set g "(JJ consists of all blocks of the following form: 



a 3 



(51) 



a n +2 1 1 



n+l 



where for 03 . . . a n +2 we take arbitrary values. If n = 0, then (ai, 02) = (1, 0), while ifn>0, ai, a2 
are determined by 



ai 



a2 



n+l ✓ \ 

= mod2 ' 

n+1 / n \ 
= 1 + ( ,-_ ! mod2 - 



(52) 
(53) 



Proofs of these propositions can be constructed using very similar argument as in Proposition 4, 
and will thus be omitted. Note that in each case, there are exactly n arbitrary values in preimages, 
so there are 2 n of them, and therefore we have 



card[g-"(? 1 )]=card[g-"(i 1 )]=2". 



(54) 



We now have enough information to obtain the number of preimages of each of the remaining 
basic blocks, since, similarly as in one dimensional case, block probabilities are interrelated via 
consistency conditions. 



Pn-l{o ) + Pn—l{ 1 ) + Pn- 



l(?0 



+ P„-l(? 1 ) = P„(0), 



P n _ 1 (0 )+P n _ 1 (%) + P n _ 1 (i ) + P n _ 1 (i ) = P n (0), 

P„-l(8 )+P„-l(8i) + i'n-l(Jo) + ^n-l(Ji) = Pn(0). 



(55) 
(56) 
(57) 



Results are summarized in Table 2, where numbers of preimages for all eight basic blocks are shown. 



b 


card[g n (b)] 


b 


card[g~ n (6)] 







2^(8-3Er= 2-^)-4-2« 


1 




2^(-4 + 4£r=o2-^)+2" 




1 




1 

1 






1 


2"V"(- 4 + 4EIU2- l( " 3) )+ 2ri 


1 

1 






1 1 


2 n 


1 

1 1 


2 "'f(4 -3 E"L„ 2" 



Table 2: Number of preimages of basic blocks for two-dimensional rule 130. 
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Dependence on the initial density in two dimensions 

Obviously, for any triangular block b we have 

p n (b)= p o(°)- 

aeg-™(6) 

Using Proposition 4 one can therefore compute P n (o 1) as follows: 



^(0l) = E E P 0(al)Pi 
i=0 a 3 ...a l+2 e{0,l} 



«3 \ 
(14 \ 

ol : jPo 

a-i+2 / 



1 

i l 



where ai,a2 are given by eq. (28) and (29). Similarly, Proposition 5 yields 



(58) 



(59) 



^(Ii) = E E M£)^o 

i=0 6 3 ...6 1+2 G{0,1} 



f>3 
6 4 



, h+2 . 



/ \ 

1 

i ..." 1 
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where 61, 62 are given by eq. (41) and (42). Now, since -P ra (l) = Pn-i{ 11) + Pn-i( )> an d 
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where, again, 01,02,61, and 62 are determined by eq. (28), (29), (41) and (42). The above is an 
exact equation of the response curve, although, unfortunately, it is not possible to calculate the 
double sums in a closed form. Nevertheless, for a given n, we can calculate and plot c n versus p, 
providing that n is not too large. This has been done for n = 12, as shown in Figure 2. Again, 
there is excellent agreement between the theoretical curve representing infinite lattice and computer 
simulations done on a finite lattice. 

It is also possible to compute P n (l) for a special case of p = 1/2. In this case Pq ( a\ ) = 1/4 for 
any a\, 02, hence 
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Figure 2: Plots of the response curve, i.e., dependence of c n on p for two-dimensional rule 130. 
Circles correspond to computer simulations, while the continuous line represents theoretical curve 
for n = 12. 



which simplifies to 
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The asymptotic density is then given by 
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where, again, it is not possible to compute the sum in a closed form. We can, however, approximate 
the infinite sum by taking the finite number of terms and use the integral bounds for the remainder, 
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The above integrals can be expressed in terms of the error function, hence for any positive integer 
k we obtain 
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where erf(x) = ^= Jq e f2 dt, and a = v/2 2 ln 2 . The sum converges very fast, and, for example, using 
the above inequalities for k = 5 we obtain 

0.1791839087 . . . < < 0.1791839597 . . . , (68) 

that is, Cqo with accuracy of seven digits after the decimal point. 
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Conclusions and future work 

We demonstrated that response curves are calculable for simple rules such as rule 130. We obtained 
exact formulae for response curves for this rule in one (eq. 20) and two dimensions (eq. 62). 
Techniques presented here are applicable to a fairly large class of CA rules - in case of elementary 
rules, we conjecture that in about 70% of cases density response could be calculated rigorously. This 
conjecture is based on numerical experiments and computerized searches for patterns in preimage 
strings. These patterns and regularities are often obvious, in other cases they can be discovered by 
constructing minimal finite state machines describing preimage sets or by using heuristic methods. 
Once the pattern is discovered, one needs to construct a formal proof that the pattern exists in 
all levels of preimages, and then perform calculations of densities. This task must obviously be 
tailored to each particular rule, and we currently do not have any general method for doing this. 
It seems, however, that rules which belong to Wolfram class I or class II are more likely to possess 
recognizable patterns in their preimage trees, and thus their density curves are often computable. 

Response curves clearly deserve further study, and it would be worthwhile to systematically 
study them for a large number of CA rules, especially in the context of applications of CA to 
solving computational problems. Work in this direction is ongoing, and will be reported elsewhere. 
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